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Abstract. 

We present an analysis of spectral line bisector variations for a few stars observed in the SARG high precision radial velocity 
planet survey, and discuss their relationship with differential radial velocities. The spectra we consider are the same used for 
determining radial velocities. The iodine cell lines employed in the measurement of radial velocities were removed before 
bisector analysis. The line bisectors were then computed from average absorption profiles obtained by cross correlation of the 
stellar spectra with a mask made from suitable lines of a solar catalog. Bisector velocity spans were then determined: errors in 
these quantities compare well with theoretical expectations based on resolution, S/N and line shape. The plot of bisector velocity 
span against radial velocity was studied to search for correlations between line asymmetries and radial velocity variations. A 
correlation was seen for HD 166435 due to stellar activity, and for HD 8071B due to spectral contamination by the companion. 
No correlation was seen for 51 Peg and p CrB, stars hosting planets. We conclude that this technique may be useful to separate 
radial velocity variations due to barycenter motion from spurious signals in spectra acquired with the iodine cell. 
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1. Introduction 

The study of activity jitter is mandatory in the search for 
exoplanets using the radial velocity (RV) technique, because it 
represents an important source of noise and a proper analysis 
is necessary to discard false alarms (e.g., HD 166435: Queloz 
et al. 2001. HD 219542B: Desidera et al. 2003 2004a_l . The 
differential RV variations induced by stellar activity are due to 
changes in the profile of spectral lines caused by the presence 
of spots and/or the alteration of the granulation pattern in 
active regions (Saar & Donahue 119971 Hatzes 120021 Saar 
1200? and Kurster et al. "2003^. The activity jitter of a star 
may be predicted by means of statistical relations from its 
chromospheric emission, rotational velocity or amplitude of 
photometric variations (Saar et al. '1998 and Paulson et al. 
|2P04i. Simultaneous determination of RV, chromospheric 
emission and/or photometry is even more powerful in disen- 
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tangling the origin of the observed RV variations Keplerian vs. 
stellar activity. However these techniques cannot be considered 
as a direct measurement of the alterations of the spectral line 
profiles that are the origin of the spurious RV variations. 
This type of study can be carried out by considering variations 
of line bisectors, that may be thought of as direct measure 
of activity jitter (Queloz et al. i2001« through the evidence of 
variations of the asymmetries which appear in spectral lines. 
Line bisectors are employed in the analysis of asymmetries 
(Gray 1982 and 1988 1 due to stellar atmospheric phenomena 
like granulation (Dravins et al. I1981> . turbulence or pulsation 
that sometimes are responsible for the RV variations observed 
in stars. 

The reason for studying simultaneously the variations of the 
RV and of the (average) line bisector of a star is to determine 
if any variation of the measured RVs are caused by a center of 
mass motion due to celestial bodies orbiting a star (see Queloz 
et al. 2001 and Desidera et al. 2003). 

The analysis of spectral line asymmetries also has an important 
application for the follow-up of transit surveys. An eclipsing 
binary in a hierarchical triple system and a transiting planetary 
companion might produce similar photometric signature; but 
very different line bisector variations, allowing one to disen- 
tangle the characteristics of the system (Torres et al. 120041 and 
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12005 > . In addition, asymmetries of spectral lines may arise due 
to contamination of the spectrum by a nearby star. This point 
is of particular relevance for the targets we are studying in the 
SARG survey (Desidera et al. 2004b ), because by design all of 
them are visual binaries. In this case a companion near the line 
of sight may contaminate the spectral features of the star being 
observed. Finally, spurious line profile asymmetries may be 
due to instrumental causes, e.g., non symmetric illumination 
of the slit. 

Line bisector variations may be studied quite easily in spectra 
acquired using fibers (see e.g., Queloz et al. l2001> . The fibers 
provide a constant, roughly symmetric illumination of the slit. 
Furthermore, spectra are generally acquired with simultaneous 
wavelength calibration lamps, rather than imprinting absorp- 
tion cell features on the stellar spectral lines. No attempt has 
been made to our knowledge to study line bisectors on spectra 
obtained through an iodine cell. One disadvantage is the 
necessity to remove the iodine lines from the stellar spectra; 
but on the other hand, the iodine lines allow a fine wavelength 
calibration and the possibility of monitoring the instrumental 
profile. 

In this paper we present such an attempt; line bisectors and 
radial velocities are determined for SARG stellar spectra, 
to study possible trends between spectral line shapes (line 
bisectors) and RVs. 

This paper describes the procedure followed to handle the 
spectra, to remove iodine lines, the construction of a mask 
made from selected lines of a solar catalogue, the cross 
correlation function (CCF) computed between mask and 
stellar spectra, the determination of the line bisector and the 
calculation of the bisector velocity span. We present results 
obtained for five stars: we found a clear correlation between 
the line bisector velocity span and RV for HD 166435, likely 
due to stellar activity, consistent with previous published work 
by Queloz et al. (2001 1. A similar correlation was found for 
HD 807 IB due to contamination by light from the companion. 
We did not find any correlation for stars known to host planets 
like 51 Peg and p CrB, and the false alarm (inconclusive) case 
ofHD219542B. 

In the next section we describe some aspects of the obser- 
vations; Section 3 explains the procedure employed in the 
analysis with results presented in Section 4. The last two 
sections are devoted to a discussion of the error analysis and 
conclusions. 
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Fig. 1. Portion of a spectrum from HD 166435. The top spec- 
trum corresponds to the program star with iodine lines, in the 
middle is the spectrum of the B-star with iodine lines and on 
the bottom is the program star spectrum after the division by 
the B-star spectrum. The division was computed for the whole 
order: the poor result is due to errors in the wavelength cali- 
bration along one order In this case the Iodine lines are better 
removed in the left half than in the right half of the order. 

of a detector with its wavelength response peaked at the 
wavelength of the iodine cell lines, was employed, keeping 
the instrumental profile stable and fairly independent of 
illumination effects. Guiding was generally done using the 
image of the binary companion on the slit viewer. 
Data reduction was performed in a standard way using IRAF' . 
High precision RVs were measured on these spectra with the 
AUSTRAL code (Endl et al. i2000: l as described in Desidera et 
al. (20031. 

The iodine cell technique includes the acquisition of spectra 
from a featureless source (a fast rotating B-star or the flat field 
lamp) taken with the iodine cell inserted in the optical path 
(Butler et al.|1996l. This kind of spectra is necessary only for 
the deconvolution of stellar templates (taken without the cell), 
but they were acquired also in most of the observing nights 
of our survey to monitor instrument performances. These 
featureless source spectra were used to remove the Iodine lines 
from the science spectra, as explained in the next section. 



2. Observations and data reduction 

The data discussed in the present paper are part of the RV 
survey aimed to find planets around stars in wide binaries 
(Desidera et al. .2004 b ). ongoing at TNG using the high 
resolution spectrograph SARG (Gratton et al. 2001 1. The 
spectra have a resolution of R ~ 150000, covering the spectral 
range 4580 A - 7900 A in 55 echelle orders, with S/N values 
in the range 70 - 400 and were taken from September 2000 
until August 2004. Slit width was set at 0.27", much smaller 
than the typical seeing Full Width Half Maximum (FWHM). 
Furthermore, an autoguider system, viewing the slit by means 



3. Data analysis 

3.1. Handling of the spectra 

The iodine lines, superimposed on the stellar spectrum for 
wavelength calibration and RV determination, were eliminated 
because the Une bisector is intended to study the asymmetries 
of the stellar spectral lines alone. For this task, spectra of fast 
rotating (Vsini > 200 km s"') B-stars were employed. The 
analysis was made only for the wavelength range where the 

' IRAF is distributed by the National Optical Observatory, which is 
operated by the Association of Universities for research in Astronomy, 
Inc., under contract with the National Science Fundation. 
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Chunk Nc. 1 



6D58 6060 6062 



6064 6066 
^length (A) 



6068 6070 



6048 6050 605? 



6058 6060 



6068 6070 5072 



6074 6076 
^length (A) 



6078 6080 



Chunk No- 6 



after application of the appropriate offset and finally divided to 
the program star spectrum (see Figure |2 details will appear 
elsewhere: Martinez Fiorenzano 120051 . The success of this 
procedure depends on the stabiUty of the instrumental profile 
over time, since the B-star spectra were usually acquired at the 
beginning or the end of each night. 

In a few cases, the cross correlation procedure used to de- 
termine the wavelength offset of chunks did not provide a 
reasonable value; in these cases the division of the star flux 
by the B-star flux added noise rather than removing the iodine 
lines. These noisy chunks were rejected from further analysis. 
Spikes due to cosmic rays or hot pixels were removed by 
replacing the spectral values within them with the averaged 
flux of adjacent spectral points. This procedure was only 
applied to those cases where spikes occurred far from the line 
centers. In the very rare cases where the spikes occurred near 
the centers of the lines used in our analysis, possibly affecting 
the line bisector analysis described below, the relevant chunks 
were simply removed from further analysis, from all spectra of 
the same star. 



Chunk Nc. 7 



6096 6098 6100 6102 6101 6106 61t 
Wo.stenglh (A) 



Fig. 2. An order of a spectrum of HD 166435 (wavelength 
range 6039 A - 6 107 A) shown in Figure [2 The order has been 
divided into 7 chunks, as described in the text. The spectra are 
shifted arbitrarily in flux to show: on top the stellar spectra with 
iodine lines, in the middle the B-star spectra (adjusted to the 
wavelength offset) with the iodine lines and on the bottom the 
clean stellar spectra after dividing the program star spectrum 
by the B-star spectrum. 



iodine lines appear in the spectra (5036 A - 6108 A) along 21 
orders. We used only this spectral region because there the 
wavelength calibration is more accurate, because of the iodine 
lines themselves. 

To handle the spectra, each order was divided in 7 pieces of 
500 spectral points each, corresponding to a wavelength width 
of ~ 10 A, overlapped by 60 points (~ 1 A) to eventually 
recover any absorption lying at the edges of the chunks. This 
procedure of cutting the spectra in pieces was clearly advan- 
tageous, while the division of complete orders displayed not 
optimal results, related to errors in the wavelength calibration 
of the spectra (see Figure 0. A cross correlation computed 
between the spectral chunks of the B-star and the program star 
determined the offset in wavelength between the iodine lines 
common to both spectra. The B-star spectrum was adjusted 
using a Hermite spline interpolation (INTEP, see Hill .1982 1 
to the wavelength scale defined by the program star spectrum 



3.2. The cross correlation function (CCF) 

3.2.1 . The solar catalogue and line selection for the 
mask 



The list of spectral lines by Moore et al. ([1966') was used to 
prepare the spectral mask needed for the computation of the 
CCF. A preliminary line list was obtained by selecting those 
lines that do not have possible contaminants (wavelength 
separation > 0.1 A), and have reduced width between 3 and 
30 F. F (Fraunhofer) is defined as the dimensionless quantity 
AA/A X lO*", where AA is the equivalent width (see Moore et al. 
1 19661 1. The range of reduced width chosen here corresponds to 
a range of central line depths from 0.14 to 0.75 in continuum 
units (in the solar spectrum). The CCFs we derived thus 
represent average profile for lines of intermediate strength. 
Once a preliminary list of lines was determined for the wave- 
length range of interest (5036 A to 6108 A), a further selection 
was made by inspecting the "The Sacramento Peak Atlas of 



the Solar Flux Spectrum" (Beckers et al. II 9761 1 . Lines were 
labeled according to their appearance: "y" (very good): sharp, 
clear, with weak wings; "y:" (good): clear but near other lines 
altering their wings; "?" (not very good): with small blends or 
strong wings. Many other lines, blended or too close to other 
spectral lines, were removed from the list at this step. 
The mask used to compute the CCF is a sum of (5-function 
profiles corresponding to 1 for the line wavelength, with a base 
of two (spectral) points, and elsewhere. The extensive line 
list (including very good, good and not very good lines: a 
total of about 300 lines) served to determine the RV; this gave 
the centroid value of the average absorption profile obtained 
by the CCF. Once this point was established, which helps to 
properly locate the profiles for the addition, a new computation 
of the CCF was performed only with the very good lines to 
obtain an improved average profile, used to calculate reliable 
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Fig. 3. An order of a spectrum of HD 166435. The left column 
shows individual chunks of the spectrum (blue lines), after re- 
moval of the iodine cell lines (see Figure Overimposed (as 
red lines) is the mask used for the determination of the CCF 
(only lines classified as very good are shown here). It is made 
of a sum of 5-functions centered on the rest wavelengths of 
the selected lines. Note the wavelength shifts between the stel- 
lar lines and the mask peaks due to the non-zero radial velocity 
of the star. The CCF for individual chunks computed by cross- 
correlating the spectra with the masks given on the left panels 
are shown in the right panels. The dashed line, shifted slightly 
below the profiles for clarity, represent the Gaussian fit com- 
puted to determine the local minimum of the CCF for each 
chunk. There are no CCF profiles for chunks 1, 4 and 7 due 
to the lack of suitable hnes for the mask in these wavelength 



average Une bisectors (see Figure [Sji. 



3.2.2. The cross correlation and addition of profiles 

A cross correlation between the mask and the stellar spectrum 
was computed for each chunk; the addition of all these cross 
correlations gave the average profile. Due to the relatively 
low S/N of some spectra, use of the average of many lines 
is appropriate for this study (the variations of line bisectors 
with time). On the other hand, the actual line bisector depends 
on the line depth, so that our "average" line bisector does not 
rigorously correspond to the line bisector of a line with similar 
depth of the CCF. Therefore, the use of our average profile 
would be misleading for some scientific goals such as the 
study of convection in stellar atmospheres. 
Due to different illumination of the CCD, each chunk along an 
order has different flux values, those close to the center of the 
orders being more luminous and yielding then results of higher 
S/N. To account for this, before summing the individual cross 
correlations, each of them was multiplied by an appropriate 
weight, proportional to the instrumental flux at the center of 
the chunk. 

Once adjusted to a common reference frame, which may be 
thought of as a centering procedure, and multiplied by its 
weight, all CCF profiles were added to get the final profile for 
the spectrum. This was normalized to determine the finally 
adopted average line profiles using the reference continuum 
determined in the IRAF reduction. This was obtained by 
interpolating a polynome throughout the spectra with a 
suitable clipping rejection procedure. Note that this reference 
continuum may contain significant errors, so that there may 
be points in the normalized CCFs that are well above unity. 
However, we kept our procedures strictly uniform throughout 
the analysis of different spectra; we then expect that these 
errors in the location of the reference continuum mainly affect 
absolute profiles and bisector estimates, but only marginally 
their spectrum-to-spectrum variations, which are of interest for 
the present discussion. 



3.3. The line bisector calculation 

The bisector of an absorption line is the middle point of the 
horizontal segment connecting points on the left and right 
sides of the profile with the same flux level. The line bisector 
is obtained by combining bisector points ranging from the core 
toward the wings of the line. 

For the bisector determination, it is necessary to adjust the 
ordinate axis of the profile to a convenient scale and to 
determine the values in velocity belonging to the left and 
right points for a given flux value. This is accomplished by 
interpolating (using INTEP, see Hill I1982> the absorption 
profile to the desired scale (details will appear in Martinez 
Fiorenzano i2005i . 
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Fig. 4. Spectrum of HD 166435. In the top panel we show the 
normaUzed cross correlation profile, the line bisector, the top 
and bottom zones (both with AF = 0.02; AF - topj - top/ = 
botf - boti). In the bottom panel we show a zoom of the profile 
with the RV scale increased to better display the asymmetries 
of the line bisector. 



3.4. The bisector velocity span 

To quantify the asymmetry of the spectral lines and look 
for correlation with RV it is useful to introduce the bisector 
velocity span (Toner & Gray il988i) . This is determined by 
considering a top zone near the wings and a bottom zone close 
to the core of the lines, which represent interesting regions to 
study the velocity given by the bisector (see Figure |4}. The 
difference of the average values of velocities in the top and 
bottom zones, Vt and Vb respectively, determine the bisector 
velocity span (BVS for short). 

The location of the top and bottom zones, as well as their width 
AF, were determined as percentages along the absorption pro- 
files. They were defined using as templates the spectra of HD 
166435; the same zones were then employed in the study of 
the other objects to be consistent in the analysis procedure. 
The Spearman nonparametric rank correlation coefficient was 
computed to establish the significance of the linear correlation 



between BVS and RV for each star. 



3.5. Error determination 



The expression given by Gray J1983I I for the bisector error due 
to photometric error was employed to establish the error for 
the bisector velocity span. The photon noise error, determined 
for the specific case of SARG using Gray's equation, yields the 
expression: 



(S\-^ AFdF\ ^'^ 

6V = { — ] \2nF 

\N \ X dV 



(1) 



where S /N is the signal-to-noise of the spectrum, including 
the contribution of the B-star spectrum and it is given by: 
((N/S)l„, + (N/S)l_^,J-^^^, with values between 65 - 300. 
The B-star spectra employed in the analysis were acquired 
on the same observing nights of the target spectra and have 
typical S/N values > 250. In eq. (1) n is the number of lines 
in the mask for the CCF; F is the central flux of the zone of 
analysis; AF is the interval of flux determining the zone; x is 
the scale of the spectrum based on the linear dispersion of the 
instrument (1.04kms"' pix"' for SARG) and dF/dV gives the 
slope of the profile in the zone being analyzed (see Figure |4ji. 
Error bars for the BVS were then determined by adding 
quadratically the errors for the top and bottom zones. 



3.6. Instrument Profile asymmetries 

Line bisector variations may arise from variations of the steflar 
spectra as well as from variations of the instrumental profile 
(IP). We studied the IP determined by AUSTRAL as part of 
the RV determination and computed line bisectors for it as for 
stellar spectra. 

Generally the line bisector variations of IP are significant, 
but they are typically much lower than the errors in the line 
bisectors determined from the stellar spectra. Corrections due 
to asymmetries of line profiles require complex convolution, 
which should take into account the real stellar profile. Actuafly, 
since the IP is much narrower than the intrinsic stellar line 
profile (the case of SARG spectra), corrections due to IPs are 
much smaller than the IP asymmetries themselves; therefore 
we did not apply corrections to the stellar BVS. 
In Figure |5]we show the run of the BVS of the IP with RV and 
the IP BVS with the stellar BVS in the analysis of HD 166435. 
No clear trend could be discerned. 



4. Results 

4.1. HD 166435 

The dwarf star HD 166435, of spectral type GO and visual mag- 
nitude my - 6.8, shows evidence of RV variations, photometric 
variability and magnetic activity. Furthermore, previous analy- 
sis of the variation of the fine bisectors revealed a correlation 
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Fig. 5. Instrumental profile (IP) analysis. Top Panel: plot of IP 
BVS vs. RVs of HD 166435. Bottom panel: plotof IPBVS vs. 
BVS of HD 166435. No evidence of trends were found. 

between RV and line bisector orientation (Queloz et al.'200ri. 
The large amplitude of the activity-induced variations and their 
stability make this star an ideal target to test new procedures 
like those presented in this paper. 

Twelve spectra of HD 166435 were acquired between May 
2003 and May 2004. The line bisectors, derived from the spec- 
tra, were used to study the correlation between the BVS and the 
RVs. The plot of BVS against RV ((Vt - Vb) vs. V,) was fitted 
to a straight line. 

In principle, the location and extension of top and bottom zones 
of the line profiles used in the derivation of the BVS are arbi- 
trary. In order to optimize detection of variations of the BVS 
correlated with RV variations, we maximized the m/cr ratio 
(where m is the slope of the straight line fitting of BVS against 
RV, and cr its uncertainty) over different choices of the loca- 
tion and amplitude of the line profile regions used to derive the 
BVS. The m/cr ratio is an indicator of the significance of the 
linear term in the straight line fitting. 

The top and bottom zones were determined according to the 
relative absorption percentages in which the highest signifi- 
cance was found: top centered at 25% of the maximum absorp- 



Fig.6. Upper panel: plot of BVS vs. RV for HD 166435 and 
the best fit to a straight line. Lower panel: line bisectors for 
individual spectra adjusted to their corresponding RV. The hor- 
izontal lines enclose the top and bottom zones considered for 
the fitting analysis. 

Table 1. Bisector velocity span from spectra of HD 166435. 53 
lines were employed in the mask for the CCF. 



JD - 2450000 


(Vt - Vb) 


Vr 


S/N 




ms"' 


m s"' 




2775.65 


-231.1 ±564 


90.6 ±5.8 


164 


2776.68 


154.1 ±57.8 


-127.9 ±6.6 


156 


2809.62 


-216.4 ±51.8 


55.1 ±8.2 


178 


2818.62 


173.8 ± 109.3 


-1044 ± 13.4 


86 


2860.40 


220.3 ± 92.8 


-83.8 ±7.7 


98 


2861.42 


110.8 ± 117.1 


32.3 ± 10.6 


81 


2862.42 


-358.8 ±64.5 


124.9 ±6.6 


140 


2891.37 


172.0 ±68.1 


-76.5 ± 7.2 


138 


2892.35 


-239.4 ± 72.5 


134.0 ±8.4 


126 


3129.74 


113.1 ±60.8 


-99.0 ± 7.4 


149 


3130.73 


76.3 ±45.9 


-16.1 ±6.1 


201 


3131.74 


-132.2 ±57.0 


70.9 ± 6.0 


164 
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Fig. 7. Upper panel: plot of BVS vs. RV for HD 8071B and the 
best fit to a straight line. The dotted line corresponds to the best 
fit discarding the spectrum with higher RV. Lower panel: line 
bisectors for individual spectra adjusted to their corresponding 
RV. The horizontal lines enclose the top and bottom zones con- 
sidered in the fitting analysis. 



Table 2. Bisector velocity span from spectra of HD 8071B. 87 
lines were employed in the mask for the CCF. 



JD - 2450000 



m s" 



Vr S/N 



1797.64 


110.3 


+ 


34.0 


-80.9 ± 


13.7 


137 


1854.54 


146.7 


± 


52.8 


28.8 ± 


11.3 


91 


2145.71 


71.6 


+ 


34.9 


-98.7 + 


11.4 


140 


2297.38 


-1.0 


+ 


44.4 


-177.7 + 


12.9 


110 


2892.61 


48.1 


+ 


42.0 


-116.6 + 


10.9 


112 


2982.48 


300.3 


+ 


62.8 


598.7 ± 


15.0 


92 


3216.73 


39.9 


+ 


33.8 


-78.7 ± 


11.0 


141 


3246.68 


82.7 
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40.3 


-74.9 


t9.2 


121 
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Fig. 8. BVS of HD 807 IB vs. the "effective contamination" 
of RV, defined as the product of the light contamination from 
the companion HD 8071 A and the RV difference between the 
components, variable with time because HD 8071 A is itself 
a spectroscopic binary. The highly discrepant spectrum from 
JD: 2452982.48 was excluded. Abscissa axis has dimensions 
of {percentage of contamination) X (m/s). 



tion, bottom at 87%; in both cases the width AF was of 25%. 
For consistency, these percentage values were then considered 
for the analysis of all the stars. 

The upper panel of Figure |6ldisplays the plot of BVS against 
RV; superimposed is the best fitting straight line, with: 

(Vt - Vb) = (-1.98 + 0.21) x V, + (-29.03 + 18.45) . (2) 

The rank correlation coefficient is -0.89 and the significance is 
of 99.99%. The lower panel of the same figure shows the line 
bisectors computed for all spectra of HD 166435, corrected for 
their RV and plotted on a common reference frame. Values of 
bisector velocity span and RV for individual spectra are listed 
in Table ^ The trend of the BVS obtained in this work is sim- 
ilar to that of Queloz et al. (2001 ); the larger value we found 
for the slope can be attributed to the higher resolution of the 
SARG spectra. 



4.2. HD 8071B 

The star HD 8071B of spectral type GOV and visual magnitude 
mv = 7.6 is a member of a visual binary system. The primary 
HD 8071A is brighter by 0.3 mag and has a projected separa- 
tion of 2.1", the smallest in the sample of the SARG survey. 
Pairs with such small separation are typically not observed in 
cases of bad seeing; furthermore the slit was oriented perpen- 
dicularly to the separation of the components. Nevertheless, 
some contamination of the spectra is still possible, especially 
for the secondary component. Our velocity measures show that 
HD 8071 A is itself a single lined spectroscopic binary with 
semiamplitude of ~ 7kms The RV difference between the 
visual components is then variable in the range (-4, 10) km s"'. 
Eight spectra of HD 807 IB were acquired between September 
2000 and August 2004. The correlation between BVS and 
RV found for all spectra has rank correlation coefficient of 
0.81 and significance of 98.5%. The analysis of the IP did not 
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show any critical trend acting on the RV computation, leading 
us to consider contamination as the cause of the observed 
correlation. The spectrum with RV ~ 600ms ' and very 
different line profile is likely heavily contaminated by the 
companion, because of the occurrence of technical problems 
related to telescope tracking that night. Our study of the line 
bisector allows us to clearly identify the problematic spectrum 
and to exclude it from the analysis of the radial velocity curve. 
Nevertheless a correlation still persists even without consid- 
ering the highly discrepant spectrum, with a rank correlation 
coefficient 0.71 and significance of 92.9%. This is likely due 
to some residual contamination (of a few percent at most) by 
the companion, compatible with the small separation of the 
pair and the actual observing conditions. Figure Q shows the 
two cases and the eight line bisectors computed for all spectra 
of HD 807 IB, corrected for their RV and plotted on a common 
reference frame. Values of bisector velocity span and RV for 
individual spectra appear in Table |2l 

To confirm our hypothesis of contamination as the source 
of observed BVS-RV correlation, we performed a simple 
modeling of the expected contamination, excluding the highly 
discrepant spectrum from JD: 2452982.48. We first determined 
the light contamination expected on the spectrum of HD 807 IB 
on the basis of the seeing conditions (given by the FWHM 
of the spectrum measured along a direction perpendicular 
to the dispersion). RV and BVS of HD 807 IB do not show 
significant correlation with such contamination. We then 
considered an "effective contamination" as the product of the 
light contamination and the variable RV difference between 
HD 8071 A and B. The Spearman rank correlation coefficient 
of RV and BVS vs. such "effective contamination" is 0.85 with 
a significance of more than 96%, supporting our hypothesis of 
light contamination as the origin of the RV and BVS variability 
of HD 807 IB. The relation between RV perturbation and 
BVS is not linear likely because the RV perturbation due 
to contamination is expected to be a non-linear function of 
the position of the contamination across the line profile (See 
Figure 
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Fig. 9. Upper panel: plot of BVS vs. RV for 5 1 Peg and the best 
fit to a straight line. Lower panel: line bisectors for individ- 
ual spectra adjusted to their corresponding RV. The horizontal 
lines enclose the top and bottom zones considered for the fitting 
analysis. 74 lines were employed in the mask for the CCF. 



4.3.51 Peg 

The star 51 Peg (HD 217014) of spectral type G2.5IVa and 
visual magnitude my - 5.5 is the first discovered solar-type 
object to host a planet, with M2 sim - 0.47 My, a = 0.052 AU 
and period of 4.23 days (Mayor & Queloz 1995 1. 
The star hes in a zone of the Hertzsprung-Russell diagram of 
stable (very low variability) objects (Eyer & Grenon 1997 1. 
Indeed, according to Henry et al. (I2000> . no measurable 
change in mean magnitude (over 5 years) was seen and the 
Ca II record displayed a signal essentially constant despite 
some season-to-season jitter and a general indication of a low 
activity level. 

Fifteen spectra of 51 Peg were acquired between June 2001 
and November 2003. There is no significant coiTelation 
between BVS and RV (rank correlation coefficient of -0.28 
and significance of 56%), confirming the results by Hatzes et 



al. ri99^) and Povich et al. ("200B. The line bisector shape 
(see Figure |5Ji seems to be constant. 



4.4. p CrB 

The star p CrB (HD 143761) of spectral type GOVa and visual 
magnitude my = 5.4 is known to host a planet with M2 sin/ 
= 1.04 My, a = 0.22 AU and a period of 39.95 days (Noyes et 
aL [T997l . 

Twenty six spectra of p CrB were acquired between April 
2001 and March 2004. In the plot of BVS against RV the 
dispersion of the data points shows no correlation (as in Povich 
et al. 1200 1> with a rank correlation coefficient of 0.15 and 
significance of 52%. 

The line bisectors and its behaviour, similar to those for 51 
Peg, are shown in Figure ^| The typical "C" shape of line 
bisectors is more evident for p CrB in agreement with the 
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Fig. 10. Upper panel: plot of BVS vs. RV for p CrB and the 

best fit to a straight line. Lower panel: line bisectors for individ- 
ual spectra adjusted to their corresponding RV. The horizontal 
lines enclose the top and bottom zones considered for the fitting 
analysis. 85 lines were employed in the mask for the CCF. 

warmer temperature of the star 



4.5. HD219542B 

The star HD 219542B is member of a wide binary system, 
with spectral type G7V and visual magnitude my - 8.6. It 
was considered in Desidera et al. (2003 1 as a candidate to 
host a planet, but ultimately discarded after further analysis; 
the small RV variations are more likely related to a moderate 
stellar activity (Desidera et al. 2004a i. 

For the present analysis of the line bisectors, only the data 
of the 2002 season were considered, twelve spectra, because 
the RV scatter and chromospheric activity were greater in 
this season (see Desidera et al. 2004a i. Therefore it should be 
easier to find a correlation between BVS and RV from these 
data alone. 

The plot of BVS against RV and the line bisectors are shown 
in Figure No correlation appears in this case: the rank 



Fig. 11. Upper panel: plot of BVS vs. RV for HB 219542B and 
the best fit to a straight line. Lower panel: line bisectors for 
individual spectra adjusted to their corresponding RV. The hor- 
izontal lines enclose the top and bottom zones considered for 
the fitting analysis. 86 lines were employed in the mask for the 
CCF 

correlation coefficient is of -0.37 and the significance of 76%. 
This lack of correlation is due to the small velocity amplitude 
(approximately between -17 and 26 ms ') and to the low S/N 
of the available spectra ~ 100. 



5. Error analysis 

An analytic study of errors was performed starting by consider- 
ing the internal BVS errors. The expected errors of BVS were 
computed for absorption profiles obtained by convolution of 
Gaussian and rotational profiles. The former were determined 
with a thermal broadening factor estimated by colors and 
temperatures of each star and the latter with a V sin i factor 
determined by the Fast Fourier Transform (FFT) analysis of 
each star's absorption profile (Grav ll992> . 
The observed errors of BVS were the quadratic mean of the 
error bars of the single spectra. The observed BVS spread was 
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Fig. 12. The expected against observed errors of BVS (tri- 
angles: quadratic mean of error bars). Squares correspond to 
observed BVS spread against the measured BVS errors with 
spread errors as error bars. 

Table 3. BVS errors: computed by eq. (1) in convolved profiles 
(expected) and measured from CCF profiles (observed). 



Star BVS error BVS error BVS spread 

observed (rms) expected observed 





ms 


ms 


ms 


-1 


HD 166435 


74.5 


85.6 


90.6" ± 


27.3 


HD 807 IB 


44.0 


48.4 


30.9" ± 


11.7 


51 Peg 


18.0 


19.1 


25.3 ± 


6.8 


p CrB 


23.0 


25.5 


28.7 ± 


5.7 


HD 219542B 


39.0 


44.0 


31.5 ± 


9.5 



These values correspond to the rms of residuals from the linear 
correlation. 



the rms average of the individual BVS for the stars that do 
not show a significant BVS vs. RV correlation (51 Peg, p CrB 
and HD 219542B), and the rms of residuals from the linear 
coiTelation for HD 166435 and HD 8071B. 
The values of errors estimated by these procedures are given in 
Table |3] The plot in Figure Elcompares the observed spread 
of data with expectations based on internal errors alone. These 
two set of values are fully consistent with each other, indicat- 
ing that our procedure is not affected by large systematic errors. 



6. Conclusions 

We studied the variation of line bisectors in the same spectra 
acquired through the iodine cell, employed also for high 
precision RV measurements. We found that such variation, 
as measured by the BVS, shows spreads fully consistent 



with internal errors, as determined from photon statistics, 
spectral resolution and intrinsic line profiles. We correlated the 
variation of the bisector span with high precision RVs for five 
stars. 

A significant correlation was established in two cases: an 
anticorrelation for HD 166435, as found by Queloz et al. 
12001 ). This is due to the stellar variability and magnetic 
activity, which makes the core of the profiles change from 
positive to negative values of RV. A positive correlation for 
HD 807 IB was due to contamination of the spectra by light 
from the companion star producing an asymmetry in the red 
wings of the profiles with a consequent inclination of the line 
bisectors toward positive values of RV. 

For the stars known to host exoplanets, 5 1 Peg and p CrB, no 
correlation was found, further supporting the conclusion that 
in these two cases RV variations are due to Keplerian motion. 
The special case of HD 219542B shows that the bisector 
technique does not allow us to disentangle the activity origin 
of low amplitude RV variations (~ 20ms"') using spectra 
of S/N ~ 100 even at the high resolution of SARG spectra 
(R ~ 150000). 

We conclude that spectra acquired using the iodine cell may be 
used to study variations of line bisectors. In order to achieve 
the required accuracy, it is necessary to deal with high quality 
spectra, in particular high S/N to reduce the error bars in the 
BVS, or to study spectra where the RV variations are of large 
amplitude. 
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